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ABSTRACT 

We report on gamma-ray observations of the supernova remnant (SNR) RX J0852. 0-4622 with the Large 
Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope. In the Fermi LAT data, we find a 
spatially extended source at the location of the SNR. The extension is consistent with the SNR size seen in other 
wavelengths such as X-rays and TeV gamma rays, leading to the identification of the gamma-ray source with 
the SNR. The spectrum is well described as a power law with a photon index of F = L85 ± 0.06 (stat);^0J9 (sys), 
which smoothly connects to the H.E.S.S. spectrum in the TeV energy band. We discuss the gamma-ray emission 
mechanism based on multiwavelength data. The broadband data can be fit well by a model in which the gamma 
rays are of hadronic origin. We also consider a scenario with inverse Compton scattering of electrons as the 
emission mechanism of the gamma rays. Although the leptonic model predicts a harder spectrum in the Fermi 
LAT energy range, the model can fit the data considering the statistical and systematic errors. 
Subject headings: acceleration of particles — gamma rays: ISM — ISM: supernova remnants — ISM: individ- 
ual objects (RX JOS 5 2. 0-4622) — radiation mechanisms: non-thermal 



1. INTRODUCTION 

Shell-type supernova remnants (SNRs) are prime can- 
didates for accelerat i on si t es of Galactic cosni i c ray s 
(iBlandford & Ostriked [19781: iMalkov & O'C Drurvl l2QQlh . 
The energy density of cosmic rays (^ 1 eV cm~^) can be ex- 
plained if ~ 10% of the kinetic energy released by supernova 
explosions is transferred to accelerated particles. Also, diffu- 
sive shock acceleration theory, which is expected to be work- 
ing at SNR shocks, can naturally explain the power-law type 
spectrum of cosmic rays. Observationally, nonthermal X-rays 
and TeV gamma rays have been detected from SNRs, partic- 
ularly from relatively young ones such as Cassiopeia A (e.g. 
iHughes et al.ll2000l:lvink & L aming 2003; Albert et al. 2007), 
RX J1713.7-3946 (e.g. > Tanaka et al. 2008: Acero et al.. 
2QQ9I: lAha ronian et al. 200 7J), and ^ J0852.0-4622 (e.g . 
Slaneetal. 2001; Ivudin etaH I2QQ5I: iPannuti et al.l 12010: 
Aharonian et al. 2007b), providing us with evidence that par- 
ticles are indeed accelerated at least up to ~ 100 TeV. 

The SNR RX J0852.0-4622 (also known as G266.2-1.2 or 
Vela Jr.) is one of the SNRs from which both nonthermal 
X-rays and TeV gamma rays are detected. The SNR was dis- 
covered in X-rays by ROSAT in the southeastern corner of 
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the Vela SNR (lAschenbachlll998h . Its X -ray spectrum was 
later co nfirmed to be of n onthermal origin by ASCA obser- 
vations ( Slane e t aH 12001 ^. The X-ray emission has an al- 
most circular morphology with a diameter of ~ 2° with rim 
brightening especially in the north, northwestern, a nd west- 
ern parts. At TeV ene rgies, the CANGAROO (Katagiri et"aL| 
2005; Eno moto et al.1 120061) and H.E.S.S. ( Aharonian eTaP 
120051 i2007bh collaborations reported the detection. Similar 
to the well-studied SNR RX J1713.7-3946, the H.E.S.S. tele- 
scope spatially resolv ed the emission from RX J0852. 0-4622 
(lAharonian et al.l[2Q0 7b) . 

The most recent and reliable estimates of age and dis- 
tance of the SNR com e from th e measurement of expan- 
sion of the shell by Katsudaetal] ^2008) using XMM-Newton 
data taken over a span of 6.5 yr. Their estimates of the 
age and distance are 1700-4300 yr and 750 pc, respectively. 
At this distance, the angular size of the SNR (~ 1° in ra- 
dius) translates to a radius of 13 pc. It is widely be- 
lieved that RX J0852. 0-4622 is a remnant of a core-collapse 
(type Ib/Ic/II) supernova explosion since studies of an X- 
ray point source, CXOU J085201. 4-461753, located close 
to th e SNR center suggest that the source is a neutron star 
(e.g. iPavlov et al.' '2 0011: iKargaltsev et al.1 l2002h . However, 
iReynoso et al. ( 2006) pointed out the X-ray source is more 
likely related to a planetary nebula they found at the location 
in radio data. 

The broadband spectrum of the nonthermal emission can 
be explained either by a leptonic scenario where the gamma 
rays are produced by inverse Compton scattering of high- 
energy electrons, or by a hadronic s cenario where tt^ decay s 
are responsible for the gamma rays (lAharonian et al.|[2007bb . 
The Large Area Telescope (LAT) aboard the Fermi Gamma- 
ray Space Telescope ( Atwood et al. 2009) can provide impor- 
tant information on GeV gamma-ray emission which is cru- 
cial to constrain the possible emission models. In this Letter, 
we report on results of gamma-ray observations of the SNR 
RX J0852.0-4622 by the Fermi LAT. 

2. OBSERVATIONS AND DATA REDUCTION 
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The Fermi Gamma-ray Space Telescope, launched on 2008 
June 11, started continuous survey-mode observations on 
2008 August 4. The LAT, with its wide field-of-view (FoV) 
of 2.4 sr, observes the whole sky in every ~ 3 hr in this 
mode. Besides the wide FoV, the LAT has better detec- 
tor performance than its predecessors, large effective area of 

8000 cm^ (on axis at 1 GeV) and point spread function 
better than 1.0° (for 68% containment) at 1 G eV. A detailed 
description of the instrument can be found in lAtwood et al.l 
(120091) . 

Here we analyzed Fermi LAT data taken for 30 months, 
from 2008 August 4 to 2011 February 8. The data analysis 
was performed using the standard analysis software package, 
the LAT Science Tools v9r21p0. We selected "Diffuse class" 
events, which are those reconstructed with high probability 
of being photons and are recommended for standard point- 
source analyses. In the analysis, the instrumental response 
functions (IRFs) "Pass6 vl 1" were used. To red uce the back- 
ground from the Earth limb (lAbdo et al.ll2009c|) . events with 
Earth zenith angle greater than 100° were eliminated. 

3. ANALYSIS AND RESULTS 

We performed a morphological and spectral analysis 
of the region that contains SNR RX J0852.0-4622. All 
the observational results except for the counts map shown 
in Figure [T] (a) and (b) are based on binned maximum 
likelihood analysis. The emission model for the maximum 
likelihood analysis includes SNR RX J0852.0-4622, Vela X, 

Eint sources listed in the First Fermi LAT (IFGL) catalog 
.bdoetal. 2010b), the Galactic diffuse background, and 
the isotropic background component which accounts for 
the extragalactic diffuse emission and residual detector 
background. One of the IFGL sources, IFGL J0854.0-4632, 
is located inside the SNR RX J0852.0-4622, and was 
excluded from the model. The spatial template for Vela X 
is assumed to be a uniform disk as reported by Abdo et al.l 
([2010a) . The Galactic diffuse and isotropic background 
models we adopted are gll_iem_v02_P6_Vl 1 _DIFFUSE.fit 
and isotropic _iem_ y02_P6_Vn_DIFFUSE. txlB, respec- 
tively. The analysis was performed on events within a 
region-of-interest (ROI) of a 14° x 14° square centered at 
(C.J2000, feooo) = (133.°00, -46.°37). 

3.1. Spatial Distribution and Identification 

In counts maps of Fermi LAT data, a gamma-ray source was 
found at the location of SNR RX J0852.0-4622. Figure[I](a) 
shows a Fermi LAT counts map of > 5 GeV events. To take 
advantage of better angular resolution and weaker Galactic 
diffuse emission, we also present a counts map of > 10 GeV 
events in Figure [T](b). No background subtraction is applied 
to the counts maps. Emission can be seen at the location of 
SNR RX J0852.0-4622. Figure[I](c) presents the test statistic 
(TS) map of the source calculated using events above 10 GeV. 
TS is defined as TS = -21n(Lo/Lpoint), where Lq and Lpoint are 
the likelihood of the null hypothesis (without a source) and 
that of a point source hypothesis, respectively (|Mattox et al.l 
[T996). To generate the TS map, we put a point source at each 
position in the sky and obtained TS for each case. The TS 
map indicates significant emission from the location coinci- 
dent with RX J0852.0-4622 above the backgrounds. 

The spatial distribution of the emission shown in Figure [T] 
appears to be extended and compatible with the SNR seen in 

http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 



Other wavelengths. To test the spatial extension and quan- 
tify it, we assumed a uniform disk as a spatial template and 
performed a ma ximum likel ihood fit to > 5 GeV events with 
pointlike (Lande "etani201 1|) . The emission was found to be 
extended significantly with TSext = -21n(Lpoint/^disk) = 144, 
where Ldisk is the likelihood for the case in which a disk is 
adopted as a spatial template for the source. The best fit ra- 
dius is l°12+^-°07 centered at (aj2ooo, W) = (133.°2 ± 0.°1, 

-U. Do 

-46.°52 ± 0.°08). The TS for the source detection is 228 with 
the best-fit disk template. The location and size are consistent 
with those of SNR RX J0852.0-4622 seen in radio. X-rays, 
and TeV gamma rays, and, therefore, support the identifica- 
tion of the emission with the S NR. We also tried fittin g the 
data using the H.E.S.S. image (lAharonian et al.ll2"007bl) as a 
spatial template. In this case, we obtained TS = 221, which is 
only slightly different from that obtained for the disk shape. 
The current statistics do not allow us to study further the mor- 
phology of the GeV emission such as spatial structures inside 
the SNR. In what follows, we use the H.E.S.S. image as the 
default spatial template to study the spectral properties of the 
source. 

3.2. Spectrum 

The spectral information of the source was extracted by per- 
forming maximum likelihood fits with gtlik^ for the ROI. In 
the energy region below 1 GeV, we found that the systematic 
errors mainly from the Galactic diffuse model become > 80% 
of the flux obtained for RX J0852.0-4622, and thus we only 
report results obtained above 1 GeV. Firstly, we analyzed all 
the events above 1 GeV together to obtain the spectral param- 
eters in the entire band. We assumed a power-law function as 
a spectral model for RX J0852. 0-4622. The IFGL sources 
except for the Vela pulsar were also modeled with power laws 
and their parameters are allowed to vary. The spectrum of 
the Vela pulsar is represented with a power law with an expo- 
nential cutoff, whose parameters (normalization, index, and 
cutoff energy) are also treated as free parameters. To the 
spectrum of the Galactic diffuse model, a power-law function 
was multiplied in order to allow for any energy-dependent 
discrepancy between the model and data. The photon in- 
dex of RX J0852.0-4622 is F = 1.85 ± 0.06 (stat)!^ }^ (sys), 
and the flux integrated between 1 GeV and 300 GeV is 
(1.01 ±0.09 (stat)!^:^^ (sys)) x 10"^ ph cm'^ s'K 

In order to obtain spectral points at each energy, we divided 
the data between 1 GeV and 300 GeV into four log-equal en- 
ergy bands, and performed maximum likelihood fits for each 
of them. The parameters for the point sources in the ROI, 
except for the Vela pulsar, are fixed at the best-fit values ob- 
tained in the entire-band fit above. The fitting parameters here 
are the normaUzation of RX J0852.0-4622, that of the Vela 
pulsar, and that of the Galactic diffuse emission. The spectral 
points obtained in each band are summarized in Table \T\ and 
are plotted in Figure [2l together with the power-law function 
from the entire band fit. 

The major contributions to the systematic errors quoted 
above come from uncertainties in the effective area cali- 
bration of the Fermi LAT and imperfect modeling of the 
Galactic diffuse emission. Uncertainties in the effective 
area are estimated based on observations of the Vela pulsar 
(^Abd o et ani2009a|) and Earth limb gamma rays (lAbdo et al.l 
2009c). Here we assumed 5% at logio(£/MeV) = 2.75, 20% 

http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
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Figure 1. (a) Counts ma p of > 5 GeV events. Th e unit of the color scale is counts per pixel with a pixel size of 0° . 1 x 0° . 1. The green contours 
are the H.E.S.S. image (lAharonian et alJlToOTbT) . Locations of IFGL catalog sources included in the emission model are indicated with the 
cyan crosses. The brightest source located to the west of the SNR is the Vela pulsar ( Ab do et al.li2009at) . The map is smoothed with a Gaussian 
kernel with a = 0?25. The inset is an image of a simulated point source with the same spectral shape as the SNR. (b) Counts map of > 10 GeV 
events. Gaussian smoothing is applied with a = 0.° 15. (c) Map of test statistic for events above 5 GeV. 



above logiQ(£^/MeV) = 4, and interpolation of the two val- 
ues in-between. The accuracy of the Galactic diffuse model 
was evaluated based on the discrepancies between the best- 
fit model and data at each location of the ROI as was done 
by I Abdo etall (l2009bD . We found that they are < 10%, and 
we changed the normalization of the Galactic diffuse model 
by ±10% from the best-fit values and added the differences 
of the parameters before and after the renormalization to the 
systematic errors. 

4. DISCUSSION 

Gamma-ray emission mechanisms for RX J0852. 0-4622 
have been actively discussed in the literature since its detec- 
tion in the TeV band. Two types of scenarios have mainly 
been considered to account for the nonthermal emission from 
RX J0852.0-4622 (e.g. Aharonian et al. 2007b). One is the 
so-called hadronic scenario in which gamma rays are domi- 
nantly radiated through decays of tt^ mesons produced in col- 
lisions between accelerated protons/nuclei and ambient gas. 
The other is the leptonic scenario in which inverse Compton 
scattering of relativistic electrons with photon fields such as 



the cosmic microwave background (CMB) is the main emis- 
sion mechanism for gamma rays. 

We calculated gamma-ray emission models for both the 
hadronic and leptonic scenarios. In the calculations, protons 
and electrons are injected with a constant rate and fixed spec- 
tral form of 



Qe,p — ^e,p ^ 



GeV c- 



exp 



(1) 



Injected particles suffer cooling via such processes as ioniza- 
tion, synchrotron, or inverse Compton scattering, and their 
spectra are deformed. The spectra of radiating particles (Ne^p) 
are obtained by solving the kinetic equations for electrons and 
protons: 



where the term be^p denotes the momentum loss rate, 
equations are solved for t = 3000 yr, and then photon spec- 
tra from electrons and protons are computed for each radi- 



(2) 



The 
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Table 1 

Fermi LAT Spectral Points 



E [GeV] dN{E)/dE [erg cm-^ s'^ ] TS 

2.04 (1.5±0.2±0.7)x 10-11 84.2 

8.49 (1.7±0.2±0.4)x 10-11 86.1 

35.3 (2.8±0.4±0.6)x 10-11 102.4 

147 (2.9 ±0.8 ±0.6) x lO-n 40.9 

Note. — The first and second errors denote sta- 
tistical and systematic errors, respectively. 



Table 2 

Parameters for the Models 



Model Sp 


PQp 


Se 


POe 


B 


n 




We 




[TeV c-i] 




[TeVc-1] 




[cm-3] 


[erg] 


[erg] 


Hadronic 1.8 


50 


1.8 


10 


100 


0.1 


5.2 X 10^0 


3.9 X 10^6 


Leptonic — 




2.15 


25 


12 


0.1 




1.1 X 10^8 



Note. — Wp and We are total kinetic energy of particles integrated above 10 MeV. 
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Figure 2. Fermi LAT spectral ene rgy distribution (SEP ) in 1- 
300 GeV with the H.E.S.S. SED by Aharonian et al.l (I2007bh plotted 
together. For the Fermi LAT points, the vertical red lines and the 
black caps represent statistical and systematic errors, respectively. 
The dotted line indicates the best-fit power law obtained from the 
maximum likelihood fit for the entire 1-300 GeV band. The butter- 
fly shape shows the 68% confidence region. The dashed curve is the 
TT^-decay spectrum by Berezhko et al. ( 2009). 

ation p rocess. The d istance to the SNR is assumed to be 
750 pc (iKatsuda et aLi r2008). As the seed photons for inverse 
Compton scattering, we considered the QMB as well as the 
interstellar radiation field in infrared and optical bands. The 
spectrum o f the interstellar radiation is taken from GALPROP 
(iPorter et al. 2006) at the location of the SNR. 

Figure [3] (a) shows the hadronic model together with ra- 
dio, X-ray, and gamma-ray data. The parameters adopted 
for the calculation are summarized in Table O The indices 
of proton and electron spectra are assumed to be the same. 
For the model curves in Figure [3] (a), the proton spectrum 
has an index of = 1.8 and a high-energy cutoff at pop = 
50 TeV c~^. In this case, the total kinetic energy of pro- 
tons above 10 MeY is Wp = 5.2 x 10^^ (n/OA cm-^y^ erg, 
where n denotes the hydrogen gas density in the post-shock 
region. The synchrotron spectrum from electrons with the 
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Figure 3. Broadband SED of RX J0852.0-4622 with (a) the hadronic 
model and (b) the leptonic model. The radio data points are inte- 
grated fluxes of the SNR determined based on the 64-m Parkes ra- 
dio telescope data by Duncan & Green (200 5). The blue line shows 
the X-ray flux estimated by Aharonian et airfcOOTbh . who analyzed 
ASCA data. The solid, dashed, and dot-dashed curves represent con- 
tributions from 7T^ decays, inverse Compton scattering, and syn- 
chrotron radiation, respectively The parameters adopted for the 
model calculation are summarized in Table [2 

same index = 1.8 is consistent with the radio spectral in- 
dex a (=-r-hl) = -0. 4 ± 0.5 de duced from the two flux 
points in Figure [3] (Du ncan & Gre en 2000). Note that the 
synchrotron spectrum index and electron index are related 
as = 2r - 1 . The flux of the synchrotron component is 
controlled by the amount of electrons and the magnetic field 
strength. The stronger the magnetic field is, the larger the 
number of electrons should be, and the inverse Compton flux 
becomes higher. In the hadronic scenario, the magnetic field 
should be at least ^50 /iG so that inverse Compton scattering 
becomes negligible in the GeV-to-TeV domain. In the case of 
B=100 fiG adopted here, the total kinetic energy of electrons 
becomes = 3.9 x 10"^^ erg, and the ratio of injected particle 
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numbers at = 1 GeV is Kgp =Ae/Ap = 1.1 x 10""^. 

In the case of the leptonic scenario, the electron index is 
constrained by the Fermi LAT spectrum. Considering inverse 
Compton scattering in the Thomson regime, the photon index 
and electron index are related as = 2r - 1 . Therefore, the ob- 
served photon index F = 1 .85 leads to a soft electron spectrum 
with Se = 2.1. With this index, however, the synchrotron spec- 
trum cannot fit the radio and X-ray data well. We then adopted 
a harder electron spectrum within the range allowed by the 
statistical and systematic errors in the Fermi LAT spectral 
points. The model curves in Figure[3](b) are the leptonic mod- 
els when the electron index is Se = 2.15 with a high-energy 
cutoff at poe = 25 TeV . The magnetic field strength can 
be determined as 5 = 12 /iG so that the synchrotron to inverse 
Compton flux ratio matches the data. Electron bremsstrahlung 
emission is calculated to be E^dN/dE ~ 10~^^ erg cm~^ s~^ 
at 10 GeV. 

Both models face difficulty. The leptonic model requires 
a weak magnetic field of order 10 /iG, which may be con- 
tradicted by observations. The width of the filaments in the 
shell observed in Chandra gives a magnetic field estimate of 
- 100 /iG (IBamba etal.1 12005: Berezhko et al.ll2009l) . These 
results are interpreted as evidence for magnetic field ampli- 
fication by streaming of accelerated cosmic rays. However, 
It is possible to reconcile a high magnetic field with the lep- 
tonic model if GeV gamma rays are radiated not only from the 
filamentary structures seen by Chandra , but also from other 
regions in the SNR where the magnetic field may be weaker. 

The hadronic model requires an unrealistically large en- 
ergy of protons if the gas density is as small as < 0.01 cm~^. 
Thermal X-ray emission works as a probe to estimate the 
gas density. However, no clear detection is reported so 
far mainly because it is difficult to separate emission from 
RX J0852.0-4622 and the Vela SNR, which overlaps and 
emits strong th ermal emission in the soft X-ray band. Us- 
ing ASCA data, ISlane et aP (1200 1|) placed an upper limit of 
n<33x 10-2 (J/750 pc)-^/^ /"^/^ cm'^ where d and / are 
the distance to RX J0852.0-4622 and the filling factor of a 
sphere taken as the emitting volume in the region used for 
their spectral analysis, respectively. In the case of / > 0.4, 
which corresponds to rather thick shell width, the upper limit 
becomes problematic since Wp reaches ^10^^ erg, which 
means that almost all the kinetic energy released by the su- 
pernova explosion must go to accelerated protons. 

The proton index = 1 .8 for the hadronic scenario is some- 
what harder than = 2.0 expected in the test-particle regime 
of diffusive shock acceleration with a strong shock. There- 
fore, it is of interest to compare the observed spectrum with 
models based on diffusive shock acceleration in the nonlin- 
ear regime. B erezhko e t al. (2009) has calculated the broad- 
band SED of RX J0852.0-4622 with their time-dependent, 
nonlinear kinetic theory for cosmic-ray acceleration. Their 
model curve is plotted over the GeV-to-TeV data in Figure O 
The Fermi LAT spectrum appears harder than the model by 
[Berezhko et al. ( 2009), but the discrepancy is not large com- 
pared to the systematic errors. 

It is of interest to compare the gamma-ray spectrum 
of RX J0852.0-4 622 with that o f the similar object, 
RX J1713.7-3946 (lAbdo et al.ll20Tl . The two SNRs have 
roughly the same age, size, and similar radio. X-ray and 
TeV gamma-ray spectra. The Fermi LAT spectrum of 
RX J17 13.7-3946 has a photon index of F = 1 .5 zb 0. 1 , which 



is more preferable for leptonic models than hadronic models 
(lAbdo etaLllfoTTh . In the case of RX J0852.0-4622, a lep- 
tonic origin is not favored but cannot be ruled out. The differ- 
ence between the two similar objects could be an important 
point to be further explored to probe gamma-ray production 
in SNRs. Further studies of Fermi LAT data in the future, 
particularly better modeling of the Galactic diffuse emission, 
will reduce the uncertainties and will allow us to probe parti- 
cle acceleration in the SNR in greater detail. 
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dation, the Swedish Research Council and the National Space 
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CNES in France for science analysis during the operations 
phase is also gratefully acknowledged. 

REFERENCES 



Abdo, A. A., et al. 2009a, ApJ, 696, 1084 

Abdo, A. A., et al. 2009b, ApJ, 706, LI 

Abdo, A. A., et al. 2009c, Phys. Rev. D, 80, 122004 

Abdo, A. A., et al. 2010a, ApJ, 713, 146 

Abdo, A. A., et al. 2010b, ApJS, 188, 405 

Abdo, A. A., et al. 2011, ApJ, 734, 28 

Acero, R, Ballet, J., Decourchelle, A., Lemoine-Goumard, M., Ortega, M., 

Giacani, E., Dubner, G., & Cassam-Chenai, G. 2009, A&A, 505, 157 
Aharonian, R, et al. 2005, A&A, 437, L7 
Aharonian, R, et al. 2007a, A&A, 464, 235 
Aharonian, R, et al. 2007b, ApJ, 661, 236 
Albert, J., et al. 2007, A&A, 474, 937 
Aschenbach, B. 1998, Nature, 396, 141 
Atwood, W. B., et al. 2009, ApJ, 697, 1071 
Bamba, A., Yamazaki, R., & Hiraga, J. S. 2005, ApJ, 632, 294 
Berezhko, E. G., Puhlhofer, G., Volk, H. J. 2009, A&A, 505, 641 
Blandford, R. D., & Ostriker, J. R 1978, ApJ, 221, L29 
Duncan, A. R., & Green, D. A. 2000, A&A, 364, 732 
Enomoto, R., et al. 2006, ApJ, 652, 1268 

Hughes, J. R, Rakowski, C. E., Burrows, D. N., & Slane, R O. 2000, ApJ, 
528, L109 

lyudin, A. R, Aschenbach, B., Becker, W., Dennerl, K., & Haberl, F. 2005, 
A&A, 429, 225 

Kargaltsev, O., Pavlov, G. G., Sanwal, D., & Garmire, G. R 2002, ApJ, 580, 
1060 

Katsuda, S., Tsunemi, H., & Mori, K. 2008, ApJ, 678, L35 
Katagiri, H., et al. 2005, ApJ, 619, L163 
Lande, J., et al. 2011, in preperation 
Malkov, M. A. 1999, ApJ, 511, L53 

Malkov, M. A., & O'C Drury, L. 2001, Reports on Progress in Physics, 64, 
429 

Mattox, J. R., et al. 1996, ApJ, 461, 396 

Pavlov, G. G., Sanwal, D., Kiziltan, B., & Garmire, G. R 2001, ApJ, 559, 
L131 

Pannuti, T. G., Allen, G. E., Filipovic, M. D., De Horta, A., Stupar, M., & 

Agrawal, R. 2010, ApJ, 721, 1492 
Porter, T. A., Moskalenko, I. V., & Strong, A. W. 2006, ApJ, 648, L29 
Reynoso, E. M., Dubner, G., Giacani, E., Johnston, S., & Green, A. J. 2006, 

A&A, 449, 243 

Slane, P., Hughes, J. P., Edgar, R. J., Plucinsky, P. P., Miyata, E., Tsunemi, 

H., & Aschenbach, B. 2001, ApJ, 548, 814 
Tanaka, T., et al. 2008, ApJ, 685, 988 
Vink, J., & Laming, J. M. 2003, ApJ, 584, 758 



